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Starting from a microscopic tight-binding model and using second order perturbation theory, we 
derive explicit expressions for the intrinsic and Rashba spin-orbit interaction induced gaps in the 
Dirac-like low-energy band structure of an isolated graphene sheet. The Rashba interaction param- 
eter is first order in the atomic carbon spin-orbit coupling strength ^ and first order in the external 
electric field E perpendicular to the graphene plane, whereas the intrinsic spin-orbit interaction 
which survives at i5 = is second order in ^. The spin-orbit terms in the low-energy effective 
Hamiltonian have the form proposed recently by Kane and Mele. Ab initio electronic structure 
calculations were performed as a partial check on the validity of the tight-binding model. 



I. INTRODUCTION 

Graphene is a two-dimensional honeycomb lattice of 
carbon atoms that has attracted considerable attention 
recently because of experimental progress 013, 0,0] that 
has raised hopes for applications in nanoelectronics and 
because of exotic chiral featuresH HItL H loL ITollTllll^ 
in its electronic structure. In the absence of spin-orbit 
interactions, the energy bands of graphene are described 
at low energies by a two-dimensional Dirac equation with 
linear dispersion centered on the hexagonal corners of the 
honeycomb lattice Brillouin zone. The recent advances in 
fabrication techniques have made it possible to produce 
graphitic systems with only a few layers or even a single 
monolayer of grapheneP,|2,l3,0|- 

One of the most remarkable properties of graphene is 
its half inte ger quantum Hall effect, confirmed by recent 
experiments [iJillJI- This electronic property follows di- 
rectly from the system's Dirac-like band structure 0, 0. 
In a recent paper, Kane and Melejj showed that sym- 
metry allowed spin-orbit interactions can generate an en- 
ergy gap and convert graphene from a two dimensional 
zero gap semiconductor to an insulator with a quantized 
spin Hall effect |1Q|. The quantized spin Hall conductiv- 
ity can be zero or nonzero, depending on the relative 
strength of intrinsic and Rashba spin-orbit interactions. 
The temperature at which the spin Hall effect can be 
observed, and the sample quality requirements for its oc- 
currence, depend on the absolute magnitude of these two 
spin-orbit interaction terms in the band structure. (Kane 
and MeleQ argued on the basis of rough estimates of the 
spin-orbit interaction scale, that the quantum spin Hall 
effect in graphene should be observable at relatively ac- 
cessible temperatures of the order of 1°K.) Motivated by 
the fundamental interest associated with the spin Hall 
effect and spin-orbit interactions in graphene, we have 
attempted to estimate, on the basis of microscopic con- 
siderations, the strength of both interactions. 

In order to allow for a Rashba interaction, we account 
for the presence of an external gate electric field E of 



the type used experimentally in graphene to move the 
Fermi energy away from the Dirac point. (Importantly 
this electric field explicitly removes inversion through the 
graphene plane from the symmetry operations of the sys- 
tem.) Then, starting from a microscopic tight-binding 
model with atomic spin-orbit interactions of strength ^, 
we use perturbation theory to derive expressions for the 
spin-orbit coupling terms that appear in the low-energy 
Hamiltonian. At leading order in ^ only the Rashba 
spin-orbit interaction term (cx E) appears. The intrinsic 
{E ~ 0) spin-orbit coupling has a leading contribution 
proportional to Both terms have the form proposed 
by Kane and Mele*^ on the basis of symmetry consider- 
ations. According to our theory the respective coupling 
constants are given by the following expressions: 
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where |s| and (spa) are tight-binding model parameters 
explained more fully below, i? is a perpendicular external 
electric field, and zq is proportional to the atomic size of 
carbon. The coupling constants Xso and have numer- 
ical values ~ 100 times smaller and ~ 100 times larger, 
respectively, than the estimates of Kane and MeleQ with 
Aso < Afi at the largest reasonable values of E. To- 
gether, these estimates suggest that the quantum spin 
Hall effect will be observable in ideal samples only at 
temperatures below ~ 0.01°K in a zero-field limit. 

Our paper is organized as follows. In Section ^] we 
briefly summarize the tight-binding model used to repre- 
sent graphene in this paper. Section ITTTI describes some 
details of the perturbation theory calculation. In Sec- 
tion IIVI we discuss ab initio density functional theory 
calculations we have carried out as a partial check on 
the tight-binding model and on the atomic approxima- 
tion for spin-orbit interactions used in the perturbation 
theory calculations. We conclude in Section with a 
brief summary and present our conclusions. 
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II. TIGHT-BINDING MODEL 
A. Two-center hopping 

For our analytic perturbation theory calculations we 
choose the simplest possible tight-binding model with 
carbon s and_p orbitals, a two-center Slater-Koster 
approximation|15j for nearest-neighbor hopping, and or- 
thogonality between Wannier functions centered on dif- 
ferent sites assumed. This gives a tight-binding Hamil- 
tonian of the form 



(3) 
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where /j,, fj! label the four orbitals on each site, A and B 
represent the two distinct sites in the honeycomb lattice 
unit cell, and Ni is one of the three vectors connecting 
a lattice site and its near neighbors. We choose a coor- 
dinate system in which the honeycomb's Bravais lattice 
has primitive vectors 



ai = a(l,0) , 
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where a = 2.46A is the lattice constant of graphene. The 
corresponding reciprocal lattice vectors are 



hi 



47r /Va 1 
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, b2 = —(0,1), 
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and the near-neighbor translation vectors are 
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The site-diagonal matrix elements are the atomic en- 
ergies of s and p orbitals, with the latter chosen as the 
zero of energy. In Table we reproduce for completeness 
the relationship between the required nearest-neighbor 
hopping matrix elements t^^^^i and the four independent 
Slater-Koster parameters [ssa), {spa), {ppa), and (pp7r) 
whose numerical values specify this model quantitatively. 
If the graphene lattice is placed in the x — y plane, Uz — 
for hops on the graphene lattice and the atomic Pz or- 
bitals decouple from other orbitals. This property is 
more general than our model, since it follows from the 
graphene plane inversion symmetry that orbitals which 
are even and odd under this symmetry operation will not 
be coupled, and is key to the way in which weak spin- 
orbit interactions influence the low-energy bands. 



B. Atomic spin-orbit interactions 

The microscopic spin-orbit interaction is 
1 
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TABLE I: Two-center matrix elements for hoping between s 
and p orbitals along a direction specified by the unit vector 
{nx,ny,n^). 



Since W is largest near the atomic nuclei, spin-orbit 
interactions are normally accurately approximated by a 
local atomic contribution of the form: 



Hso — E ■ 



(8) 



where j is a site index, Pu denotes projection onto an- 
gular momentum I on site i, ^; is the atomic spin-orbit 
coupling constant for angular momentum I, and S is the 
spin operator on site i. For our model spin-orbit coupling 
occurs only among the p orbitals. 



C. External gate electric fields 

Finite carrier densities have been generated in 
graphene by applying an external gate voltage. The re- 
sulting electric field E lifts inversion symmetry in the 
graphene plane. An electric field E can also be pro- 
duced by accidental doping in the substrate or cap layer 
or by atomic length scale charge rearrangements near the 
graphene/substrate or graphene/cap-layer interfaces. To 
model this important effect we consider an additional lo- 
cal atomic single-particle Stark-effect term of the form 



H 



EF 



^eEY^: 
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where i is a site index. In our s — p tight-binding model 
the only nonvanishing matrix element of Hef is the one 
between the s and Pz orbitals to which we assign the 
value cEzq. 



III. PERTURBATION THEORY 

A. Unperturbed Hamiltonian matrix at K and K' 

The low-energy Hamiltonian is specified by the Dirac 
Hamiltonian and by the spin-orbit coupling terms at K 
and K' . We choose the inequivalent hexagonal corner 
wave vectors K and K' to he K = i(26i + h) = (|f , 0) 
and K' ^ -K. Table HD and Table iHll list the Hamil- 
tonian matrix elements and the corresponding eigenvec- 
tors. Here s is the on-site energy of s orbitals relative 
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to p orbitals, a = ^{spa), /3 = | [{ppa) — {ppir)], and 
_ Vs'^+8a^±s ^ Note that the a bands are decoupled 
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TABLE II: Tight-binding model matrix elements at the K 
and K' points in the absence of spin-orbit interactions and 
external electric fields. The first (second) sign corresponds to 
the K{K') point. 



from the tt bands. When the spin-degree of freedom is 
included, the E = Q eigenstates at K and K' are fourfold 
degenerate. Below we refer to this degenerate manifold 
as D. 
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TABLE III: Unnormalized unperturbed eigenvectors at the K 
and K' points arranged in increasing order of energies assum- 
ing < 7+ < 2/3 < 7_. The first (second) sign corresponds 
to the K{K') point. 
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TABLE IV: The effective spin-orbit matrix at the K point. 



Similar results are obtained at the K' point. It follows 
that the effective spin-orbit interaction for tt orbitals is 

HefS — ~y^SO + ^SO CTzTzSz + {(TxTzSy — ffySx) (12) 

where the CTq. Pauli matrices act in the A, B space with az 
eigenstates localized on a definite site, — ±1 for K, K' 
points, and the Sa are Pauli matrices acting on the elec- 
tron's spin. This Hamiltonian differs from the form pro- 
posed by Kane and MeleQ only by the constant —\so- 
The excitation spectrum has a gap Egap = 2(Aso — ^r) 
and the system has a quantized spin Hall effect |^ for 
< Ai^ < Aso- 

To obtain quantitative estimates for the coupling con- 
stants we used the tight-binding model parameters listed 
in Table IVI taken from Ref. For the spin-orbit cou- 
pling parameter among the p orbitals we use ^ = 6 meV, 
a value obtained by fitting carbon atomic energy levels 
given by the ah initio electronic structure code described 
below. These values imply a graphene energy gap at 
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B. Low-energy effective Hamiltonian 

We treat the atomic spin-orbit interaction and the ex- 
ternal electric fields as a perturbation: 



AH — Hso + Hef- 



(10) 



The effective Hamiltonian which lifts the E — degen- 
eracy is given by the second-order degenerate state per- 
turbation theory expression j^^: 
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where m,n Cz D. An elementary calculation then shows 
that the matrix elements of Hm'n (at the K point) are 
those listed in Table lTvl with A50 and \r defined by Eqs. 
and ((2Jl respectively. 



TABLE V: Hopping parameters for a graphene taken from 
Ref. m 



Ai^ = equal to 



2A5C 
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« 0.00114 meV « fcs x 0.0132°K, 



(13) 

where we used the nonorthogonal tight-binding parame- 
ters neglecting the overlap for simple estimations. Our 
estimates of \r are discussed later. 



IV. AB INITIO CALCULATIONS 

We have performed realistic ah initio electronic struc- 
ture calculations ^3 inversion symmetric (A^; — 
0) graphene sheets using the projector augmented 
wave (PAW) |3 method with a Perdew-Burke-Ernzerhof 
(PBE) generalized gradient approximation (GGA)(l9l| 
density functional in order to partly test the quantitative 
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accuracy of the conclusions reached here about spin-orbit 
interaction gaps based on a simpHfied electronic structure 
model. The calculations were performed using VASP (Vi- 
enna ab initio simulation package) 20j. In VASP, spin- 
orbit interactions are implemented in the PAW method 
which is based on a transformation that maps all elec- 
tron wave functions to smooth pseudo wave functions. 
All physical properties are evaluated using pseudo wave 
functions. The spin-orbit interaction is evaluated tak- 
ing into account only the spherical part of the potential 
inside mufhn tins surrounding the carbon nuclei: 



demonstrate that the tight-binding model spin-orbit gap 
estimates are accurate. 
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In order to make the gaps induced by spin-orbit interac- 
tion exceed the accuracy of VASP eigenvalues, we have 
artificially increased the strength of Hso by up to 300 
times by decreasing the speed of light c. 

Figure ^ shows the tight-binding band structure of 
graphene for ^ = and ^ = 300^0, where = 6 meV. 
The spin-orbit gap is not large on the scale of the full 
band width, even when enlarged by a factor of 300. 
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FIG. 1: (Color online) Graphene band structure for ^ = 
and f — 300^0 using the tight-binding model with nonorthog- 
onal orbitals. Hopping parameters were taken from Ref. 2V 
and ^0 = 6 meV was used for the atomic spin-orbit coupling 
strength. 

Figure |21 compares the ab initio calculation and 
tight-binding model low-energy gaps at the hexagonal 
Brillouin-zone corners for Xj^ = 0, finding close agree- 
ment. Both approximations find a gap that grows as 
the second power of the spin-orbit coupling strength. 
The close agreement is perhaps not surprising given that 
VASP also makes an atomiclike approximation for the 
spin-orbit coupling strength. In our opinion, however, 
the neglected contributions from interstitial regions and 
from aspherical potentials inside the muffin-tin sphere are 
small and their contributions to energy levels tends to- 
ward even smaller values due to spatial averaging by the 
Bloch wave functions. We believe that these calculations 
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FIG. 2: (Color online) Energy gap for = as a func- 
tion of spin-orbit coupling strength from the ab initio cal- 
culation, from the tight-binding model with nonorthogonal 
orbitals, and from the analytic expression in Eq.Q. 



V. DISCUSSION AND SUMMARY 

The intrinsic and Rashba spin-orbit interactions arise 
from mixing between vr and a bands due to atomic spin- 
orbit interactions alone in the case of Xso (Eq-CQ) and 
due to a combination of atomic spin-orbit and Stark 
interactions in the case of Xr (Eq.Q). These expres- 
sions for Xso and follow directly from Ea. ((TTl) and 
from the eigenvectors and eigenenergies listed in Table 
mil (The energetic ordering in Table IIIII applies for 
< 7+ < 2/3 < 7_ which holds for the tight-binding pa- 
rameters in Table Ivl) The pure p — p hybridized bonding 
and antibonding states (energies ±2/3 in Table llllj) are 
symmetrically spaced with respect to the undoped Fermi 
level and do not make a net contribution to either X^i or 
Xso- The s —p hybridized bonding states (energy — 7_ 
in Table ITTl|l . on the other hand, are further from the 
Fermi energy than the corresponding antibonding states 
(energy -I-7+ in Table Illlj) because of the difference be- 
tween atomic s and p energies. Their net contribution to 
Aso is proportional to s and inversely related to (spa), 
which sets the scale of the energy denominators. Sim- 
ilar considerations explain the expression for Xr which 
is proportional to ^ and cEzq and inversely proportional 
to spa. Xr vanishes at £^ = because of the inversion 
symmetry of an isolated graphene plane. 

The numerical value of the Rashba interaction param- 
eter Xr obviously depends on the electric field perpen- 
dicular to the graphene plane which varies as the carrier 
density is modulated by a gate voltage. A typical value 
can be crudely estimated from Eq. by assuming a 
typical electric field E 50V/300nm5], and using the 
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value zq Sas x (O.620I/O.529A) obtained by scaling 
the hydrogenic orbital Stark matrix element by the ratio 
of the atomic radii |23 of carbon and hydrogen: 

Xfi = g ^ 0.0111 meV fcg x 0.129°K. (15) 

6[sp(T) 

The value of Xj^ is influenced by screening of the elec- 
tric field at one graphene atom by the polarization of 
other graphene atoms and by dielectric screening in the 
substrate and cap layers, but these correction factors are 
expected to be '--^ 1. Note that our estimate for Xso is 100 
times smaller than Kane's estimate, ~ 1°K, whereas Xn 
is 100 times larger than Kane's estimate, ~ 0.5°mK. If 
our estimates are accurate, Xso < Xr at large gate volt- 
ages. For undoped samples, however, the requirement 
for a quantized spin Hall effect gap0, Xso > A_r, should 
still be achievable if accidental doping in the substrate 
and cap layer can be limited. When Xso is smaller than 
Xr, the energy gap closes and graphene becomes a zero 
gap semiconductor with quadratic dispersionjSi]. 

Our estimates suggest that the quantum spin Hall ef- 
fect in graphene should occur only below ~ 0.01°K, a 
temperature that is still accessible experimentally but 
not as convenient as ~ 1°K. In addition, it seems likely 
that disorder will dominate over the spin-orbit couplings 
in current samples, so further progress in increasing the 
mobility of graphene sheets may also be necessary be- 
fore the quantum spin Hall effect can be realized experi- 



mentally. We emphasize, however^^, that the spin Hall 
effect survives, albeit with a reduced magnitude, even 
when the spin-orbit gap is closed by disorder. 

In summary, we have derived analytic expressions for 
the intrinsic and Rashba spin-orbit interaction coupling 
constants that appear in the low-energy Hamiltonian of 
a graphene sheet under a perpendicular external electric 
field. The Rashba interaction parameter is first order 
in the atomic carbon spin-orbit coupling strength ^ and 
the perpendicular external electric field E, whereas the 
intrinsic spin-orbit interaction is second order in ^ and 
independent of E. The estimated energy gap for E — 
is of the order of 0.01°K and agrees with realistic ab initio 
electronic structure calculations. 

Note added in proof. Recently, we became aware of 
two other articles which address spin-orbit interactions 
in graphene and reach broadly similar conclusions [2^. 
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